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Abstract. We construct C°° solutions to the one-dimensional nonlinear wave 
equation 

utt — Uxx — -^^i— ^ = with p > 
that blow up on any prescribed uniformly space-like C°° hypersurface. As a 
corollary, we show that smooth solutions can blow up (at the first instant) on 
an arbitrary compact set. 

We also construct solutions that blow up on general space-like C'° hyper- 
surfaces, but only when 4/p is not an integer and k > (3p -|- 4)/p. 



1. Introduction 

In this paper we consider the one-dimensional focusing nonhnear wave equation 

Uu-u,.-^^\u\Pu^O (1.1) 

for general powers p > 0. The coefficient of the nonlinearity can be changed via 
the rescaling ux{t,x) := u{Xt^Xx). However, this particular choice simplifies many 
of the formulae that follow. 

Even for smooth initial data, solutions to (|l.ll) can blow up in finite time. The 
simplest example of this is 

u{t, x) = {T - ty'^/P defined for t<T and x eR. (1.2) 

(Note the benefit of choosing the coefficient 2{p + 2) / in (|1.1|) .) 

The equation (II. ip is relativistically invariant. This blurs the distinction between 
space and time, as well as implying finite speed of propagation. As a result, it is 
natural to consider the maximal space-time region on which the solution is defined. 
Moreover, this region must take the form 

{{t,x) : a{x) < t < a{x)}, 

where either a(x) = cx) or ct is a 1-Lipschitz function, that is, 

sup '"(f-"f"^' <l, 

x,y \x-y\ 

and either cr(x) = — oo or fj is a 1-Lipschitz function. (The number 1 represents 
the speed of light in our model.) A point t = cr{x) on the (forward) blowup surface 
is called non- characteristic if 

\a(x) — a(y)\ 
limsup J— V < 1; 

v^x \x-y\ 

otherwise, it is called characteristic. The surface {t = cr{x)} is called space-like if 
it is everywhere non-characteristic. 

1 



2 



ROWAN KILLIP AND MONICA VI§AN 



In this paper we discuss the question of which surfaces can occur as blowup 
surfaces. In the papers Oil], CafFarelli and Friedman showed that a must be at 
non-characteristic points for a closely related model (they also obtain conditional 
results in higher dimensions). More recently, Merle and Zaag considered precisely 
the model (|l.ip in a series of papers; see [lOl [11] and the references therein. They 
have shown that the characteristic points form a discrete set and that a is away 
from these points. (The also obtain a wealth of information about how the solution 
behaves near the blowup surface.) In [T^], Nouaili extended their methods to show 
that away from the characteristic points, a is C^'" for some (unspecified) Holder 
exponent a > 0. 

In the other direction, one may seek to show that certain rich classes of functions 
a do indeed arise as blowup surfaces. The only works we know of that address 
this type of question are the papers [U [9] of Kichenassamy and Liftman. These 
authors consider arbitrary spatial dimension with p e N and show that any analytic 
space-like surface can occur as a blowup surface. In his lecture notes on the subject, 
Alinhac explicitly asks about the possibility of extending the work of Kichenassamy 
and Liftman to the C°° case; see [21 P- 9]. In this paper, we work only in one 
dimension and show that general C°° space-like surfaces can indeed occur as blowup 
surfaces. Our methods apply also in the case of finite regularity; however, for 
technical reasons, the argument becomes messier if 4/p is an integer, so we are 
omitting that case. 

Theorem 1.1 (Regular solutions that blow up on a space-like hypersurface). 
Suppose either (i) p > and k = oo, or (ii) ^ ^ N and 3-|-| <A:€NU {oo}. Let 
'E — {t = cr{x)} be a uniformly space-like hypersurface [i.e., sup^, |cr'(a;)| < 1). 
Then there exist a nearby space-like hypersurface Eq := {t = cro(a;)} with a{){x) < 
a{x) and a solution u G C'l^^{{ao{x) < t < o-{x)}) to (jl.l|) . with m — [fc — 3 — 4/pJ , 
that blows up on S in the sense that 



for all X eR. 

Remarks. 1. The floor notation, [a:J , denotes the greatest integer less than or equal 
to X. Similarly, the ceiling \x~\ is the least integer greater than or equal to x. 

2. The blowup rate in ()1.3|) is Lorentz invariant: to this degree of accuracy, 
u{t, x) agrees with the Lorentz boost of the solution (|1.2p that blows up on the 
tangent line to E passing through the point t — a{x). 

3. The space-like hypersurface {t = ao{x)} may follow the blowup surface {t = 
cr(a;)} so closely that our solution u is not defined globally in space at any one 
particular time. (This issue also appears in results of 9 .) 

4. The condition fc > 3 -I- 1 arises as a natural requirement for our method. This 
is discussed at some length in Section \5\ 

Our initial interest in treating the case of smooth (as opposed to analytic) blowup 
surfaces was to address the question alluded to in the title: Can solutions to nonlin- 
ear wave equations blow up on Cantor-like sets? By this we mean that at the time 
when the first singularity occurs, can this singularity arise simultaneously across a 
Cantor-like set? Note that if a is an analytic function, then {x : cr{x) = T} must be 
discrete (or all of M); however, given an arbitrary compact set E, it is not difficult 




(1.3) 
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to construct a smooth function for which E is a. level set. In this way, Theorem 1 1.1 1 
allows us to construct a solution to (|l.ip that blows up on an arbitrary compact set. 
Moreover, we are able to control the time of existence sufficiently well that we can 
guarantee that the solution admits global (in space) Cauchy data (cf. Remark 3 
above). The precise statement is as follows: 

Corollary 1.2 (Blowup on arbitrary compact sets). Fix p > 0. Given a compact 
set E G there exist C°°(R) initial data (u(0),ut(0)) and a smooth solution to 
(jl.ll) with this initial data that blows up at a time < T < oo precisely on the 
set E, which is to say, 

• u{t, x) converges to a smooth function of x £ R \ E when t ^ T and 

• livcit^T u(t, x) = oo for all x € E. 

The proof of Theorem I 1 . 1 [ proceeds by constructing a parametrix ii for the desired 
solution and then using contraction mapping arguments to correct it to an actual 
solution u to (jl.ip . Due to the singularity of the solution we are seeking, the 
equation obeyed by the correction term u — £t is of Fuchsian type, which is a PDE 
analogue of an ODE having a regular singular point. This has two consequences. 
First, the parametrix needs to be constructed to very high order. This is done in 
Section [3l Second, the contraction mapping argument revolves around a singular 
variant of the linear wave equation. The basic estimates for this equation are 
derived in Section |3] and then used to construct the solution in Section [5] 

Working in one spatial dimension gives rise to a number of simplifications, both 
mathematical and expository. We make particular use of the Minkowski-space ana- 
logue of the Riemann mapping theorem (Proposition 12.21 below) . This allows us 
to choose coordinates {s,y) for the space-time region {t < a{x)} that straighten 
out the boundary while only deforming the d'Alembertian, — 9^, by a conformal 
factor. That this is only possible in R^+^ is a famous discovery of Liouville; specifi- 
cally, the only (Euclidean or Minkowski) conformal mappings in higher dimensions 
are Mobius transformations (see O §15]). As the notion of a conformal mapping 
relative to a pseudo-Riemannian metric is discussed in relatively few differential 
geometry textbooks ([5 is one exception), we develop the facts we require from 
scratch in Section [21 

While Mobius transformations in Minkowski space have played a central role 
in several fundamental developments in the theory of nonlinear wave equations 
(spectacularly so at the hands of C. Morawetz or D. Christodoulou, for example), 
we have not seen the greater generality available in the M.^'^^ setting used previously 
to aid in the analysis of this type of equation. 

Acknowledgements. We are grateful to Mario Bonk for valuable discussions 
about conformal mappings in Minkowski space. 

The first author was partially supported by NSF grant DMS-1001531. The 
second author was partially supported by the Sloan Foundation and NSF grant 
DMS-0901166. This work was completed while the second author was a Harrington 
Faculty Fellow at the University of Texas at Austin. 

2. Existence and properties of conformal maps 

In this section, we use a conformal mapping (relative to the Minkowski metric) to 
reduce the problem of constructing solutions to (jl.ip that blow up on a (sufficiently 
smooth) space-like hypersurface to that of constructing solutions that blow up on 
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the spatial axis; these latter solutions must solve a suitably modified nonlinear wave 
equation. 

Definition 2.1 (Conformal mapping). A diffeomorphism $ : {s,y) i"> {t,x) 
between two open sets in Minkowski space R^+^ is called conformal if there exists 
a conformal factor A(s, y) > such that 



dt 

ds 



ds 



at 

dy 
dx 
dy. 







§1 

ds 

M. 

dy 



ds 
dx 
dy. 



= A(s,y) 



-1 




(2.1) 



Equivalently, the pull back of the Minkowski metric via $ differs from the original 
Minkowski metric by the scalar factor A(s, y) > 0. The mapping $ is a (Minkowski) 
isometry at those points where A(s, y) = 1. 

From (I2.ip it is clear that the null (or light-like) directions are unchanged by a 
conformal mapping. As A > 0, the notions of space- and time-like are preserved. 
Thus there are four classes of conformal mappings depending on whether the spatial 
and/or temporal orientations are preserved or inverted. When both are preserved, 
the mapping is called proper (cf. [5] §6.2]). 

It is easy to derive the Minkowski analogue of the Cauchy-Riemann equations 
from (|2.1I) . In the case of proper mappings, they read 



and # = 1^. 

ay OS 



In particular, (t, x) are wave coordinates, that is, tgs — tyy — and Xss ~ Xyy ~ 
(in the sense of distributions). Using the general solution of the one-dimensional 
wave equation, we deduce that all proper conformal mappings of Minkowski space 
take the form 



t=W{y + s)~g{y-s)] and x ^ ^[fiy + s) + g{y ^ s)] , 



(2.2) 



where / and g are orientation-preserving diffeomorphisms of R. (The re- 
maining three types of conformal mappings correspond to allowing / and/or g 
to be orientation-reversing.) A more geometrical derivation of this fact (which 
we learned from Mario Bonk) is as follows: As conformal mappings preserve the 
null directions, the grid of light rays in (i, x) coordinates must be mapped onto 
the corresponding grid in {s,y). Thus, a conformal mapping corresponds to an 
(arbitrary and independent) change of variables in each of the null directions: 
(x + t, X ^ t) = {f{y + s), giy — s)), which is, of course, just another way of writing 

(El. 



Proposition 2.2 (Existence of conformal maps). Let 2 < k E NUjoo} and consider 
a uniformly space-like hypersurface S = {t = (j{xY\. Then there exists a 
conformal diffeomorphism from {s > 0, y G M} onto {t < a{x)}, which is an 
isometry between the boundaries. Moreover, the conformal factor A(s,y) is C^y^. 

Proof. The natural conformal mapping $ : (s, y) i— ;> [t, x) is not proper. We seek 
instead one for which (s,y) i— )■ $(— s,y) is proper, which, from the preceding dis- 
cussion, must take the form 



x + t = f{y — s) and x — t — g{y + s) 
with increasing functions f,gEC'^ (R) . 



(2.3) 
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As E is uniformly space-like, |c'(x)| < 1 uniformly for x G M. Thus, by the 
implicit function theorem, there exists a C*^ function : M — > K so that 

t — a{x) if and only if x + t ~ h{x — t). 

Moreover, < h'{z) = < e'^ for some e > 0. 

The condition that the conformal mapping sends {s = 0, y G R} onto S reads 

f(y) = h{g{y)) for ah y e M, 
while the condition that it is an isometry on the boundary is 

ny)9'{y) - 1 for ah y G R. 
Therefore, g must solve the following ODE: 

9'{y)^[h'{9{y))Y^'^ with, say, .g(0) = 0. 

(We take the positive square-root.) By Picard's existence and uniqueness theorem, 
this equation admits a unique global solution g ^ for h G with k > 2. Thus 
f ^ h o g £ C'^ , as well. Moreover g'{y) > e and /'(y) > £, so we are guaranteed 
that both are diffeomorphisms onto the whole of R. 

A computation shows that the conformal factor is given explicitly by 

Xis,y) :^f'{y-s)g'iy + s)>0, 

which is clearly C^y^. It is now an easy matter to verify that the conformal mapping 
we constructed satisfies the desired properties. This completes the proof of the 
proposition. □ 



A direct consequence of Proposition 12.21 is the following 

Corollary 2.3 (Reduction to blowup on lines). For any 2 < fc G NU{oo}, there is a 
one-to-one correspondence between solutions to (jl.ll) that blowup on a C'^ uniformly 
space-like hypersurface in the sense of (jl.3l) and solutions to 



Vss - vyy - ^^A(s, y)\v\Pv = (2.4) 
that blow up on {s = 0, y G M} in the sense that 

\im spv{s,y) — 1 (2-5) 

s\0 



for each y G R. Here, A(s,y) denotes the conformal factor from Proposition[ 

Proof. Let <i> : (s,y) i— > it,x) be the conformal mapping from Proposition! 
straightforward computation shows that 

d',^dl^Xis,y){d?~d',), 

which means that v{s,y) — no $(s, y) solves (12. 4p if and only if u{t, x) solves 

To prove the correspondence between the blowup behavior of v on the line {s = 
0, y G R} and that of u on the hypersurface S, we merely have to show that 

lim '^^''^ ~ * = -Vl - \<j'{x)\^ = lim "^^^^ ~ ^ (2.6) 

i/'cr(x) S s\0 S 

where {t^x) = $(s,y). In the former limit, x is held fixed; in the latter, y is fixed. 
To verify these relations, we need to compute the derivative of $ (and its inverse) 
at the boundary. 
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Recall from the proof of Proposition 12.21 that our conformal mapping takes the 
form 

t = ^[fiy - s) - g{y + s)] and x = ^[f{y - s) + g{y 
with / and g increasing and obeying 

1 r. , ..-1 1/2 



f'iy) 



9'{y) 



[h'{g{y))\ 



-l+V(£)]l/2 
.l-cy'{x)\ 



when X — \ [f{y) + g{y)\^ that is, when s = 0. Therefore, 



-dt 


dt_- 


ds 


9y 


dx 


dx 


.ds 


dy. 



-f'{y)-9'{y) f{y)-9'{y) 
-r{y)+9'{y) f'iy)+9'iy) 



1 



-1 a'ix) 
^a'ix) 1 



where, as before, (0,t/) and (cr(x),x) = $(0,7/) are corresponding points on the 
boundaries. The left column of this matrix justifies the second identity in (|2.6|) . 
while the top left entry of its inverse explains the first identity. Note that this 
matrix has determinant —1, as can be predicted from the knowledge that A = 1 at 
the boundary and $ is orientation-reversing in the first variable and orientation- 
preserving in the second. □ 

3. Building the parametrix 

In this section, we further analyze the modified wave equation 

Vss-Vyy-^^Xis,y)\v\Pv^O, (3.1) 

where A denotes the conformal factor introduced in the previous section. We remind 
the reader that the requirement that the conformal mapping be an isometry on the 
boundary translates to X{0,y) = 1 for all y G M. Our goal is to find a good 
parametrix for (|3.ip that blows up on the line {s — 0, y G K}. 

Let us start by recalling a very simple version of the well-known extension the- 
orems of Whitney and Borel (cf. Theorem 1.2.6 and Corollary 1.3.4 in |6 ): 

Lemma 3.1 (Whitney/Borel extension). Given k G {0, 1, . . . , cx)} and functions 
Pj{y) G Cy^^M.) for < j < k + 1, there exists a function p G C^y(R'^) obeying 

[dlp\{S),y) = p,{y) for each 0<j<fc + l. 
Naturally, if k = oo, then k — j = k + I = oo. 

We will use Lemma 13.11 to construct the parametrix when ^ is not an integer. 
The remaining case is slightly more involved and will be addressed afterwards. 

Proposition 3.2 (Parametrix). Let p > such that | is not an integer. Fix 
2< k eNU {oo} and A e C^-^{[0,oo) x M) such that X{0,y) = 1 for all y £ R. 
Then there exists a nowhere-vanishing parametrix v G C^y^ {{0,oo) x M) for the 
equation p.ip whose failure to satisfy the equation can be described as follows: 
setting 

£ := Vss - Vyy - ^^|t^A(s, y)\v\''v, 

we have 

}+k-3-a\ 



dfdP£{s,y)^o{s- 



0, 



(3.2) 
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for all y (z R and for allO < a + /3 < k. When k = oo, this means that the left-hand 
side of p.2p vanishes to all orders as s — > 0. Moreover, v blows up on the line 
{s = 0, 2/ G M} in the sense that 

2 

\im spv{s,y) = 1 for all y£R. (3-3) 

Proof. We start by noting that if A = 1, then there is a weU-known solution to (|3.1I) 
that blows up on the line {s = 0, y G IR.}, namely, v{s,y) = s"^/''. Therefore, we 
will look for a parametrix of the form v(s,y) = s~'^^Pp{s,y) with p{0,y) = 1. The 
condition that v satisfies p.ll) to the level of accuracy described in p. 21) translates 
to the following condition for p: 

d:d^£{s,y)^o{s'^-'-"), (3.4) 

uniformly in y G M and for all < a + /3 < fc, where 

S ■■= Pss - Pyy - Is-'ps + ^S'^ip - X\p\Pp]. (3.5) 

This in turn corresponds to conditions on Pj{y) [9J/o](0, y) for < j < fc. 

We first seek to satisfy p.4p in the case a — (3 — 0; that (I3.4p continues to hold 
for all < a + /3 < fc will be an immediate consequence of our construction. As 
p(0, y) = 1, we must have po{y) = 1, which ensures cancellation of the 0{s~^) term 
on the left-hand side of (j3.4p when a — 13 = 0. To continue, we expand A in Taylor 
series near s to maximal order, 

fc— 1 

X{s, y) = 1 + sdAiO, y) + • • • + j^zrjy di'-'^m v) + o{s^-'), 

and we substitute this into (I3.5P . With a little computation, one finds the following 
coefficients for on the left-hand side of p.4p (in the case a = /? = 0): 

ifj = -2 



hiiv) - ^dsX[Q.y) ifj = -l 



.HMf)p^(y)_HM, 

I (ItS - |)^.+2(y) - j^dlp.iy) + P, if J > 0, 

where Pj is a polynomial in dsX{Q,y), . . . , 9^"'"^A(0, y) and pi{y), ■ ■ ■ ,pj+i{y). 
To zero the coefficient of s~^, we simply choose 

Piiy) -f±|9.A(0, y)eC^^ 

More generally, since 4/p is not an integer, the factor in front of pj^2{y) is non- 
zero and so we may recursively choose this function to annihilate the 0{s^) term. 
Simple induction shows that pj e C^~^~^ . The continuation of this procedure is 
limited only by the regularity of A. More precisely, when k is finite, the process 
terminates with the choice of pk-i G C*'(R) to eliminate the 0{s^~^) term on the 
right-hand side of p. 51) . 

By Lemma [STTl there exists p G C5^^^([0,oo) x R) which matches our Taylor 
expansion at s = 0. This guarantees that p.4[) holds, which in turn shows that the 
parametrix v G C^_^i((0,oo) x K) satisfies and ((H^)) . 

To see that the parametrix v can be chosen to be nowhere-vanishing, we note 
that as poiy) = 1, the extension p{s,y) can be chosen to be nowhere-vanishing, 
indeed, bounded from below by a positive value. This can been seen quite directly 
from the proof of Lemma 13.11 or by post facto modification of any other extension. 

This completes the proof of the proposition. □ 
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The proof above breaks down in the case 4/p is an integer: it is impossible to 
cancel the ©(s"*/^) term in E because the coefficient of P2+i/p vanishes. This phe- 
nomenon is well-known in the context of power-series solutions to linear ODEs and 
our treatment of this case is adapted from that model. For simplicity of exposition, 
we will only consider the case where A is Cf'y. 

To construct the parametrix in the case when ^ is an integer (and the conformal 
factor is smooth), we will rely on the following analogue of the Borel extension 
lemma, whose proof is easily adapted from that of [BJ Theorem 1.2.6] cited earlier: 

Lemma 3.3 (Borel extension). Given C^(M) functions Pj/{y) for integers < 
j < oo and < i < L{j) < oo, there is a function p G C^((0, oo) x R) obeying 



V ^ 

j=0 £=0 



for any < a, /3, J < oo and any e > 0. 

Proposition 3.4 (Parametrix). Suppose | G N and A G Cf^([0,oo) x R) with 
A(0,2/) = 1. Then there is a nowhere-vanishing parametrix v G C^((0,oo) x R) for 
the equation ()3.ip such that (I3.3|l holds and 

£ := Vss - Vyy ^ A(s, ?/) | w j^u 

along with all of its derivatives vanish to infinite order as s 0. 

Proof. Consider what happens when we substitute 

1+^ 



pToiy) , P2+^,liy) 2+1, . ^ 



3 = 



^ jl (2 + ^ 

^ p 



into p.5|) . To make the coefficients of s^, for — 2 < £ < in the resulting expansion 
vanish, we choose pj,o to be the functions pj from Proposition 13.21 for all < j < 
1 -I- ^. Moreover, the coefficient of s'*/^ log(s) vanishes for the same reason that 
the coefficient of P2+4/p vanished in the proof of Proposition 13.21 It is possible to 
choose P2+4/p,i to cancel the coefficient of s**/^, because this coefficient is now 

(3 + ^)P2+l,i(2/) d^p.^^oiy) 



(2 + |)! (|)! 



P4 



where, just as in the preceding proof, P^/p is a polynomial in pi(y), . . . , Pi+4/piy), 
as well as the first 2 + A/p derivatives of A. 

To continue, we recursively choose coefficients in the expansion of the type given 
in Lemma l3.3l to construct our full parametrix. The success of this procedure relies 
on the fact that for j > A/p, one can choose Pj+2.1 to cancel the term of order 
s-' log^(s) in the expansion of £, where £ is as in (|3.5p . To see that this is indeed 
possible, we note that the coefficient at this order is 

where Pj^i is a polynomial in the Taylor coefficients of A and lower-order coefficients 
in the expansion of p. By lower-order coefficients we mean pji^t with j' < j or j — j' 
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and i' < £. Note also that for our application, L{j) — for < j < I + ^ and more 
generally, L{j) = 0{j) for any < j < oo. 

That p(s, y), and hence the parametrix w, can be chosen to be nowhere- vanishing 
follows from the fact that po,o(y) = 1 and L{0) = 0. This completes the proof of 
the proposition. □ 



4. Basic estimates for a singular wave equation 

In this section, we prove finite speed of propagation and basic estimates for 
solutions to the critically singular wave equation 

Vss - Vyy - {v^ - \)s-'^v = F. (4.1) 

In our applications, v ~ 3/2 + 2/p. The reason behind the appearance of the 
singular perturbation will become apparent in Section [SJ see (j5.3p . 



Proposition 4.1 (Estimates for (|4.ip ). Fix v > and let k{s,y;so) be the causal 
fundamental solution to (|4.1I) . that is, 

kss - kyy - {v"^ - \)s^'^k = 5{s - SQ)5{y) for s, sq G (0, cxi) and y e R (4.2) 

with fc(s, y; sq) = for s < sq. Then k obeys finite speed of propagation, namely, 

k{s,y;so) = when \y\>s — so, (4-3) 

and satisfies 

||fc(s,2/;so)|Li(K) < (s-so)(^)'' (4.4) 
Moreover, dsgk and dgk are finite measures with total variation 

\\dsok{s,y;so)\\ + \\dsk{s,y; so)\\ <1 + — ", (4.5) 



and, when v 



\\dsk{s,y-s^)\\^,^ <l+'-^(lX'\ (4.6) 

When V ~ 1/2, k is the fundamental solution for the regular wave equation, 
namely, 

, I ^ 1 1 f c \ sin(2<[s - So]) 
fco(s, y] So) = -jX{\y\<s-sa} Or fco(s, so) = • 

The properties (|4.3p . (|4.4p . and (|4.5p are obvious in this case. 

Just as for the usual wave equation, the fundamental solution k carries all the 
information necessary for the Duhamel formula: 

v{s) ^ k{s;so) *Vs{sq) - ■§^{s;so) *v{so) + k{s; si) * F{si) dsi, (4.7) 



whenever < sq < s and v solves (j4.ip . Here * represents convolution in the 
y variable. The estimates appearing in Proposition 14.11 allow us to address the 
possibility of sending sq i in this formula, which we will do in the proof of 
Corollary 1121 
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Proof. Fix > ^ . Claim (14.31) is inherited from the usual wave equation, as can be 
seen by iterated Duhamel expansion. We will also give an analytical proof of this 
in what follows. 

As equation (|4.1I) is invariant under spatial translation, one can compute k via 
its spatial Fourier transform 

k{s,^-so) := / e-2"^«fc(s,y;so)rf2/. 



Taking the spatial Fourier transform in (14.21) we obtain 

kss + 47r2^2fc - {u^ - i)s"2fc = 5{s - sq). 
This is a version of the Bessel equation and has solution (cf. 1 , §9.1.49]) 

_ /fV^[J.(27r|C|so)y.(2^|^|5)-r.(2^|e|so)J.(2^|ek)], s > 



k{s,C,so) 



0, s < So, 



with Ji, and Y^, denoting the Bessel functions of the first and second kind, respec- 
tively. Notice that fc is a continuous function in s, while dgk has a jump discontinuity 
at s = So! indeed, dsk ^ I as s sq, as follows from the Wronskian relation for 
Bessel functions (cf. [U §9.1.16]). 

For fixed s, sq £ (0, oo), k extends from ^ > to an entire function in C which 
is even under the transformation ^ i— > — ^. Indeed, while both J^, and Y^, typically 
have branch singularities at the origin, these cancel out in the combination that 
forms k; see [D §9.1.2, 9.1.10, 9.1.11]. 

One can also expand k in terms of the Hankel functions hI^^ = J,y + iY^ and 
Hi/ ~ J u lYi/'. 

u c , _ /|iV^[^fi'^(27r^so)-ffi'^(27res)-Hi'^(27r^so)i^i''(2^es)], s > sq 
[0, s < So. 

Using the behavior of the Hankel functions in the whole complex plane (cf. [TJ 
§9.2.3, 9.2.4]), we find that 

|Ms,^;so)hO.,.o(e'"l''"^ll'"^''') as |C| ^ oo. 

Thus, by deforming the contour in the imaginary direction (up if y > 0, down if 
2/ < 0), one can deduce 

k{s,y;so)= f e^^'y^k{s,^;so)d^ = for \y\ > s - so > 0. 

JR 



This completes the proof of 

It remains to prove (|4.4p . (|4.5p . and (|4.6I) . As the kernel kg associated with the 
usual wave equation obeys the bound (j4.4p . we need only control on the difference 
k — ko in Ly. Both kernels also obey (|4.3I) and so 



^~^o\\rUK} - ~ *o|||fc - fcoll ^oo|-R^ < 2|s - solllfc - fc| 



ri, 



This is the strategy we use to prove (|4.4p ; compare (|4.12p below. We will prove (|4.5p 
in an analogous manner; however, a more careful subtraction needs to be made. 
Specifically, we will show that 

||a.„fc(s,f; so) - a,„fco(s,C; so) + ^^kois,^-, so)\\^,^^^ < So"'(^)'"^- (4-8) 
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The analogue of (|4.8p for dsk is 

||9,Ms,e;so)-9.fco(s,C;so)-^^fco(s,e;so)||^i(j,) <s-i(^)''"\ (4.9) 

We will not give details of the proof of (|4.9I) , because it differs very little from that 
of (|4.8p . Note that (|4.6p follows from (14.91) : the necessity of taking v > ^ comes 
from estimating the contribution of fcg. 

To complete the proof of (j4.4p and (14. 5p . we need only crude bounds on the 
Bessel functions for small values of ^, namely, 

|J^(z)| < min(2'',z-i/2j < max(z-^z-l/2) (410) 

for each fixed v > ^ and uniformly in z £ (0, 00). For larger f we use 

Mz)+^Y4z)=Hi'Hz) ^ {^f\,iz)e^^-^-^''^+'^/\ (4.11) 

where z S [l,oo) and is complex-valued and obeys (cf. [1, §9.2.5-10]) 

a^{z) = 1 + + 0{z-^) = 1 + 0{z-'^) uniformly for z G [I, 00). 

Remember that and are real- valued for such z. 

To control dsgk, we combine the above with the relations 

d.Mz) = f J,(z) - J,+i(z) and d,Y,{z) = i^Y,{z) - r,+i(z); 

see [H §9.1.27]. This yields 

a,„fc(s,e;so) = fVii^{^[J.(2^|e|so)>^(2^|e|5) -r.(27r|C|so)J.(27r|e|s)] 

-27r|C|[J.+i(27r|e|so)y.(27r|e|s)-y,+i(2^|e|.so)J.(2^|^|s)]}, 

whenever s > sq. 

We now turn to the details of showing (|4.8I) and 

||Ms,e;so)-fco(s,e;so)||^i(R) < (^)' (4.12) 

As some terms can be treated simultaneously, we use a G {0, 1} to indicate the 
number of sq derivatives. We divide the integrals into three regions, as follows: 
Region I: |^| < . For ^ in this region, a G {0, 1}, and s > sq, we have 

|5?o^(S''?;so)| + |<9",fco(s,^;so)| < — Sg", 

as follows easily from (j4.10p . Thus, we can estimate the contributions of this region 

by 

|9«fc(s,e;so)|-f |af„fco(s,e;so)|dC< (4.13) 
Region II: 3-^ < \^\ < Sg ^ Using (HTU)) again, we find 

Therefore, we bound the contribution of this region by 

/ \df^kis,^;so)\ + \d^M',^;so)\dC<s-"(-y'\ (4.14) 
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In the trivial case i/ — ^, integrating the bound given above actuaUy produces a 
factor of log(s/so) rather than (s/sq)" = 1. 

Region III: |^| > Sq ^. Now we use formula (|4.1ip . This quickly yields 

|Ms>C;so) - fco(s,C;so)| < —TjT^ 

and consequently, 

|fc(s,C;so)-fco(s,f;so)|de<l< (-) \ (4.15) 
For the derivative, the analogue of the first calculation is longer; it yields 
\dsMs,^;so) - a.„fco(s,C;so) + ^^fco(s,e;so)| < 

(The O(ICr^) terms collapse so elegantly because '^^"^^ll'"^ " ^ = ^^fe^') ^his 
leads immediately to 



l«l>So' 



e; so) - 5«„fco(5, so) + ^^fco(s,C; so)| < So \ (4.16) 



which is consistent with (|4.8p . 

Combining (|iT^ . (|iTil) . (|iT5|) . and (|i?TBl) proves and (|i?T^ . which in 

turn complete the proof of the proposition. □ 

Corollary 4.2 (Estimates on solutions to the singular wave equation). Fix v > ^, 

6 > v—^, and a compact interval G / C [0, oo). Consider the initial-value problem 

\vss - Vyy - (i^^ - j)s^'^v = F 
\^;(0)=0, v,{0)^0 

for V : I X R ^ R. If F e C^yil x M) with s~^F{s,y) e Lf^y{I x M), then this 
admits a distributional solution with v and Vg continuous; moreover, for integers 
a, /3 > and s G /, 

fj.=0 

whenever the right-hand side is finite. 

Proof. If F''^^ G C^y and supported on [e, oo) x M, then the standard theory of 
strictly hyperbolic equations applies (cf. |7j Ch. XXIII]) and yields a solution 
v'^^'> G C^y that vanishes when s < e and also obeys 

[d's - - (-^ - i)s-^]afa>(^) = d^Fi^^ 

Q-l 

where Fjf^ dfF^^^ + {v^ - i) ^ (^) [^r'^s"'] [^^^^^j ■ 

Ai=0 

The key issue is to obtain bounds independent of e and convergence as e ^ without 
relying on the additional smoothness of F'^^^ 
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For 0<e<s<e'<oo, 



< 



k{s, z; so)d!^F^'> {sa,y- z) dz dso 



and so by (14.41) we see that for (5 > — |, 

a-l 

s.V s,y 

Combining this with induction in a yields 



uniformly for 0<e<s<e'<c» and ?/ e M. In addition to showing convergence, 
these bounds also provide uniform control by sending e J, and e' — > oo. 

By using (j4.6p in addition to the above, it is possible to control one more deriv- 
ative of v^'^^ : 



df+'d^v'^'Hs,y)\ < 



ks{s, z; so)d!^Fjf\so,y - z) dz ds[ 



Notice that the additional derivative costs one power of s on the right-hand side. 

It remains to estimate u^^^ — v^^ ^ in the region 0<e<£'<s< oo. By the 
Duhamel formula (|4.7p . 



+ r\kis;so) * d^M'H-^o) - Fpiso)]\dso. 



Combining (|4.4p and (|4.5p . with 

\d:d^v(^'\e')\ <F (e' - ef{e''f and \df+^d^v^'\e')\ <f (e' - e){e')\ 

which where proved previously, we see that the first two terms obey 

\k{s; e') * ar^a.^^^' (£') I + I ^ (s; e') * {e') I 5"^+^ [e'f-''+He' ~ e), 

and hence converge to zero as e' — >■ 0, since S > v — ^. 

To estimate the contribution of the inhomogeneity, we argue exactly as in the 
case £ < s < e' to obtain 

f]k{s;so) * d^[Fi-\so)-Fi-'\so)]\dso < s^+'J2P^'''''9^^^^ ' 

where ly — ^ < S < S and Lfy relates to the slab [e', s] x R. Using the fact that 
S < S when s is small and that i^^^' — to cover larger values of s, we see that 
the L°° norm above converges to zero as e' — 0. 

This completes the proof of the convergence and thus settles the lemma. □ 
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5. Construction of the blowup solution 

In this section we prove Theorem 11.11 and Corollary 1 1.2 1 Theorem 11.11 will follow 
easily from Corollary 12. 31 and the following 

Theorem 5.1 (Regular solutions to ()3.1|) that blow up on a line). Suppose either (i) 
p>0 andk = oo, or (ii) | ^ N anrf 3 + 1 < fc e NU{cx)}. Let X G C^-^{[0,oo)xR) 
such that A(0, y) — 1 for all j/ e M. Then there exist a time So{p,X) > and a 

solution V € ci'^y''^^^^^^ ((0, sq] x ^) to p.ip which blows up on the line {s — 0, y e 
R} in the sense that 

lim spi;(s, y) = 1 for all y € R. (5.1) 

s\,0 

Additionally, Sq admits a positive lower bound that depends only on the cl,^^^^^^ 
norm of A. 

Proof. Given sq, (5 > and an integer m > 0, consider the Banach space X™{[0, sq]) 
consisting of those w E C^y{[0, so] x R) for which the norm 

is finite. This is the family of spaces in which we will run a Picard iteration 
(— contraction mapping) argument to improve the parametrices v from Proposi- 
tions 13.21 and 13.41 to an actual solution v — v + w. In terms of the spaces just 
introduced, the error in the parametrices, namely, 

£ Vss ~ Vyy - ^^^^^A(s, y)\v\Pv, 

satisfies 

ll'^ll-^r-i([o,'*o]) - 1 whenever m + 5 < k-3~ ^. (5.2) 

In fact, this quantity is o(l) as Sq — > 0; however, we will not make use of this. 
In order for v = v + w to he a solution to p.l[) . we must have 

Wss - Wyy - A(s, y){\v\Pv - \v\Pv) + £ = 0. 

By construction, v{s,y) = s~^^Pp{s,y) with p € C^^^ and p(0,y) = 1. Using this 
and Taylor's formula, the equation for w reduces to 

Wss - Wyy - (i^^ - ij'S^^ic = F, (5.3) 

where v — 1 + -, 

Fis, y) := -i + 'J£±^£±lls-'{XpP ~ l)w + ^J£^X f" {i + Ow) (1 - 6) dO, 

and f{w) — \w\Pw. We will soon estimate F in the space Xg^^. The motivation 
for this (and indeed, the introduction of these norms themselves) is the fact that 
Corollarv 14.21 m'ovides the necessary estimates for solutions to (|5.3p . Specifically, 

\\w\\xpi[o,so]) Sn so\\F\\x^_^([a^so]) provided 5 > |. (5.4) 

Note that to derive this, one must consider (|4.17p with varying values for S, not 
merely that appearing in the formula above. 
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To continue, note that {XpP){0, y) — 1 and so, with a httlc extra work, one finds 
\\d^d^s-\XpP <s-^~^' when fi + l3<k-l. 

Therefore, 

\\s-^iXpP - j[o,so]) ^ ll^ll^rdo^'^o]) provided m < fc - 1. (5.5) 

We now move to the third term in F, which involves f{w) — Iwl^w. Note that 
/ fails to be C°° at w ^ 0. This will not be a problem for us, because for sq small 
we have v{s, y) ~ s~'^^p uniformly on [0, so] x moreover, from the argument that 
follows, we see that ||iy||L~j^([o,so]xR) ^1- As a result of this, 

|/(^+2)(^) + ^y;)|_ s-l(p-^-i) on [0,so]xM 

uniformly for 6 E [0, 1]. Using this we find 

Q+/3 I 

i=0 aha 0\-f^ j=l 

where cJ, /3 are partitions of a, /3 into £ parts. This directly yields 

d"dyl"{v + Ow) < fi-f (P-i)-" exp(||w||x-) for a + /? < m and any S > 0. 

Naturally, the implicit constant depends on m. Proceeding from here we find 

\\Xw^f"{i + 0w) 11;^™ ^([o,,„]) ^ \Mx',-([oM) exp(||u)||x™([o,so])) (5-6) 

uniformly for 9 g [0,1], provided m < k. Indeed, one can exhibit additional 
(helpful) positive powers of sq on the right-hand side. 

With minor modifications, these estimates also show that 

< \\w - w\\x^iioM)e:ip{2\\w\\x^Hio,so]) +'2\\w\\x^^(lo,so]))- 
Let us now put all the pieces together. We choose m< [fc — 3 — |J, finite when 
fc = oo, and (5 > | so that m + 6 < fc — 3— |. When fc < oo this is possible, since in 
this case A/p is not an integer. By combining (|5.2p . (|5.4p . (|5.5p . and (|5.6p . we see 
that the sequence of Picard iterates w„ obey 

\\Wn+l\\x^([OM) - •5o{l + ||w„||xj-([o,.o]) + ll"'nllx,"([0,.o])^^P(ll^"ll^r([0,^o]))}- 

Thus by choosing so sufficiently small, we can guarantee that Wn are uniformly 
smah in ([0,so]). 

The term estimated in (15.51) is linear in w and so the bound works equally well 
for differences. Using this and (|5.7p we obtain 

||w„+i - w„||x^" < so||w„ - exp(2||ti;„||x^r„ + 2||w„_i||xjr.), 

which shows that the Picard iterates do indeed converge (for sq small). 

For each integer m < [fc — 3— -J, we have proved the existence of a C™ 
solution to p.ip on a small neighborhood of s = 0, which additionally satisfies 
(|5.ip . Naively, the interval [0, sq] of existence may depend on m. However, simple 
(inductive) persistence of regularity arguments using only Duhamel's formula for 
the one-dimensional wave equation show that the solution v cannot blow up in 
C^y while remaining bounded in C^ y. We use, of course, that A € C'^,y^- This 
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observation has two important consequences: First, in the case fc = oo it shows 
that the solution we construct is indeed C°°. Second, it shows that the time of 
existence depends only on the cl"^^^^^^ norm of A. □ 

The reader will have undoubtedly noticed that in order to produce any solution 
at all (say merely C^^y), this argument requires a great deal of regularity on the 
conformal factor A. This is not purely an artifact of the proof, but reflects a 
fundamental underlying phenomenon, which we explain next. 

Even in the case A = 1, there are multiple solutions that diverge at s = at 
exactly the rate prescribed in (|5.ip . Indeed, there is a one-parameter family of such 
solutions that depend on s alone. To be more precise, using energy conservation, 
one sees that any such solution v{s) obeys 

l + 2^iJa-(P+2)]-^da, 

v(s) 

where E = — ^jt'l^^^ denotes the (conserved) energy density. From this it 

is easy to derive the asymptotic behavior 

Thus, to uniquely specify such a solution we require a parametrix v for which 



which corresponds to a = and k = [3 + |] in (|3.2p . In this sense, our regu- 
larity requirements are minimal for the construction of w by contraction mapping 
arguments, indeed, by any method that ensures uniqueness. 

We now return to the main theme of the section and prove the results stated in 
the introduction. 



Proof of Theorem As mentioned at the beginning of the section, the construc- 
tion will be a direct application of Corollary 12.31 and Theorem 15.11 Suppose there- 
fore that either (i) p > and A; = cx), or (ii) | ^ N and 3-|-|<fceNU {oo}. 
Let T, = {t — cr(a;)} be a C'^ uniformly space-like hypersurface. Then by Corol- 
lary 12.31 the solution v to (j3.1|) described in Theorem 15.11 gives rise to a solu- 

[fe— 3— - J 

tion u e Ct ^ " to (dH]) in a neighborhood of the hypersurface, namely, on 
{tTo(a^) < ^ < ""(x)}, where the Cauchy hypersurface Yjq = {t — (Tq{x)} is the 
image of the line {s = sa(j), X), y G M} via the conformal mapping $ from Proposi- 
tion 12.21 Here so(p, A) is as in Theorem 15. II Moreover, by Corollarv 12.31 and ()5.ip . 
the solution u has the desired blowup behavior. □ 



Proof of Corollarv \1.2[ Let C M be a compact set. We start by constructing 
a smooth function that has E a.s a, level set. Indeed, let Ij denote the bounded 
connected components of R \ i? and let : R — > [0,1] be a smooth function, with 
(j){x) = when a; < 0, 4>{x) > when x > 0, and 4>{x) = 1 when x > 1/2. We 
define the hypersurface T, = {t = a{x)} as follows: 

X eE 

aix) ^e{l + exp(-^)0(iiHM)), x e I, 

(/)(dist(a;,£')), X < inf E or X > sup E. 
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Here e > is a small constant that will be chosen later. Observe that ||cr||c'= ^fc £j 
for each integer k. 

Now consider the conformal mapping $ of {s > 0,y G R} onto {t < <j{x)} 
constructed in Proposition [2]2] From this construction, we see that for e small, 

||^*-Id|L. <ke 

for any integer k, where denotes the derivative of $. Thus, 

dist($"i({0,x}),{s = 0}) < £ 

and the conformal factor obeys 

IIA-lllc- <ke 

for each integer k. Combining these two observations, we see that when s is suffi- 
ciently small, the smooth solution to 

V.S - Vyy - ^^X\V\PV = 

that was constructed in Theorem 15.11 will exist throughout the region $^^({0 < 
t < a{x)}). (Of course, this relies on the fact that so(p, A) depends robustly on A.) 
Therefore, u = wo$~^ is a smooth solution to defined on {0 < t < (j{x)} that 
blows up at time t = e precisely on the set i?. □ 
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